Epidemics in small world networks 

M. M. Telo da Gama and A. Nunes 
Centra de Fisica Teorica e Computacional and Departamento de Fisica, 
Faculdade de Ciencias da Universidade de Lishoa, 
P- 1649-003 Lisboa Codex, Portugal 

Abstract 

For many infectious diseases, a small- world network on an underlying regular lattice is a suitable 
simplified model for the contact structure of the host population. It is well known that the contact 
network, described in this setting by a single parameter, the small-world parameter p, plays an 
important role both in the short term and in the long term dynamics of epidemic spread. We have 
studied the effect of the network structure on models of immune for life diseases and found that 
in addition to the reduction of the effective transmission rate, through the screening of infectives, 
spatial correlations may strongly enhance the stochastic fluctuations. As a consequence, time series 
of unforced Susceptible-Exposed-Infected-Recovered (SEIR) models provide patterns of recurrent 
epidemics with realistic amplitudes, suggesting that these models together with complex networks 
of contacts are the key ingredients to describe the prevaccination dynamical patterns of diseases 
such as measles and pertussis. We have also studied the role of the host contact strucuture in 
pathogen antigenic variation, through its effect on the final outcome of an invasion by a viral strain 
of a population where a very similar virus is endemic. Similar viral strains are modelled by the 
same infection and reinfection parameters, and by a given degree of cross immunity that represents 
the antigenic distance between the competing strains. We have found, somewhat surprisingly, that 
clustering on the network decreases the potential to sustain pathogen diversity. 

PACS numbers: 89.75.-k, 87.10. +c, 87.23.-n 
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I. INTRODUCTION 



Nearly a decade ago, Watts and Strogatz introduced a class of networks with a topology 
interpolating between that of lattices and random graphs P] . In these models a fraction of 
the links of the lattice is randomized by connecting nodes, with probability p, with any other 
node. For a range of p the network exhibits 'small world' behaviour, where a local neigh- 
bourhood (as in lattices) coexists with a short average path length (as in random graphs). 
Analysis of real networks reveals the existence of small worlds in many interaction net- 
works, including networks of social contacts. Recently, attention has been focussed on the 
impact of network topology on the dynamics of the processes running on it with emphasis on 
the spreading of infectious diseases. Earlvnumerical and analytical studies have focussed on 
the calculation of epidemic thresholds 14 M for different classes of networks and on the short 

Q 

term dynamics of epidemic bursts j6l| for small world networks of the Watts and Strogatz 
type. By contrast, the effects of network topology on the long term dynamics of epidemic 
spread has received much less attention. 

In the last couple of years we have studied simple epidemic models on dynamic (or 
annealed) Watts and Strogatz small world networks [7J and found that, for a large range 
of p, epidemics persist for times that approach those of homogeneously mixed populations 
providing suitable models to investigate the effects of spatial correlations on the long term 
dynamics of epidemic spread. In our model the N nodes of the network represent individuals. 
The links are then connections along which the infection spreads. If all the links are local 
lattice) the persistence of the epidemic is low and the endemic state may not be reached 
a]. On the other hand, if the network is random (homogeneously mixed population), the 
persistence is high but the effect of spatial correlations is neglected ^. The dynamic version 
of the Watts and Strogatz model [7 1 is thus the simplest model of a real network of social 
contacts 

Using a SEIR (Susceptible-Exposed-Infective-Recovered) model, we have obtained a 
quantitative description of the incidence oscillations that characterize the prevaccination 
recurrent epidemics of immune for life diseases such as measles and pertussis. The model, 
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on the same network, was found to describe the quahtative differences observed in the long 
term dynamics of such diseases. We have also used a simple two-strain version of a recently 
proposed SIR model with reinfection to determine the outcome of an invasion in multi- 
ple strain infections. The results uncovered a surprising effect of the network structure on 
the diversity of pathogens, that was found to decrease, rather than increase, as the spatial 
correlations on the dynamic small world network increase. 

II. RECURRENT EPIDEMICS ON NETWORKS 

The Susceptible-Infected-Recovered (SIR) mean field model in a closed population of 
individuals is given by the equations 

S = -pSI/N 

i = f3SI/N--fI (1) 

where S (respectively /, R) is the number of susceptible (respectively infected, recovered) 
individuals and (3 is the transmissibility. This is a model for the dynamics of a disease that 
confers permanent immunity and it is easily generalized to an epidemic that takes place 
on a network. As for models of epidemic spread on regular lattices jsl, the model may be 
mapped onto percolation on the same network j^. The percolation transition corresponds 
to the epidemic threshold, above which an epidemic outbreak is possible (i.e. one that 
infects a non-zero fraction of the population, in the limit of large populations) and the size 
of the percolating cluster above this transition corresponds to the size of the epidemic. The 
SEIR model is a generalisation of that takes into account a latency period, during which 
individuals carry the pathogen but are not yet infective. This was found to be relevant in 
the modelling of childhood diseases j^. 

The endemic persistence of a disease on a closed population requires renewal of suscepti- 
bles through demography or loss of immunity. To assess the effect of network topology on the 
long term dynamics of simple spatial models for childhood diseases we implemented discrete 
SIR and SEIR models with births and deaths corresponding to an average life time of 61 
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lOl lll | on a square lattice with sites and dynamic small world interaction rules. 
We found characteristic medium and long-term dynamics related, in a quantitative fashion, 
to the structure of the network of contacts. In particular, as the small world parameter, 
p, decreases, the increase in spatial correlations (i) decreases the effective transmissibility 
through the screening of infectives and susceptibles, which in turn increases the value of 
the transmissibility at the endemic threshold. In addition, the spatial correlations (ii) en- 
hance the stochastic fluctuations with respect to the homogeneously mixed stochastic model. 
This effect is particularly strong at low p, where the relative fluctuations are largest and 
where as a consequence (iii) the dependence of the steady state densities on the effective 
transmissibility predicted by the mean- field equations breaks down j^. 

In Figure 1(a) we have plotted the incidence time series of a SEIR model with the 
epidemiological parameters (latency period ti and recovery time Tj) taken for measles and 
two different values of the small world parameter p, p = 1 and p = 0.2. For homogeneously 
mixed populations, p = 1, the infective time series exhibits incidence oscillations with an 
average period of two years. However, the amplitude of these oscillations is underestimated 
when compared with the incidence oscillations of measles in real data for similar population 
sizes ^2^. By contrast, in the infective time series for SEIR simulations on the network for 
p = 0.2 the amplitude of the incidence peaks is shown to increase significantly, in line with 
the real data. 

As with any unforced model, the characteristic biennial cycles of prevaccination records 
cannot be obtained without fine tuning of the model's parameters, transmissibility, (3, and 
probability of long-range infection, p. Despite this fact, which means that, for measles, 
seasonality cannot be ignored the reported effects of spatial correlations should also be 
taken into account as a key issue in explaining incidence amplitudes. 

Figure 1(b) shows the incidence time series of model with the epidemiological param- 
eters taken for pertussis and three different values of the small world parameter p, p = 1, 
p = 0.2 and p = 0.1. Again, we see that for the smaller value of p the amplitude of the inci- 
dence peaks is significantly enhanced with respect to the homogeneously mixed case p = 1. 
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However, for intermediate values of p, and in particular for p = 0.2, this effect is almost 
negligible. 

The main conclusion of the previous analysis is that the impact of the spatial correlations 
depends also on the stability of a particular disease, measured by the distance of the endemic 
phase to the corresponding threshold on the network, that is disease dependent jQ]. 



III. EVOLUTIONARY DYNAMICS OF MULTIPLE STRAINS 

We consider the invasion by a pathogen strain of a population where another strain is 
endemic. A single viral strain is meant to represent a cluster of co-circulating antigenically 
close strains, and its uncoupled dynamics is described by a SIR model with transmissibility 
f3 for first infection and partial immunity against reinfection. The reinfection parameter a of 
the model is the ratio between disease transmissibilities for reinfection and for first infection. 
It measures the size of the cluster, in terms of antigenic distance. 

For realistic values of the birth rate, the SIR model with partial immunity exhibits an 
abrupt change in the levels of infection as f3 increases across a transmissibility threshold. 
This region of abrupt change of infection prevalence is related to a bifurcation that occurs 
in the model with zero birth rate |7], and we shall use the term reinfection threshold for 
ihe smoothed transition that takes place in models with small birth rates. Recent studies 
2U W — a between the .eintectton t..e=Md and the potential tot 

pathogen diversity, that increases as the level of infection increases. 

In two strain dynamics, disease transmissibility for infection by one of the strains of 
individuals previously infected by the other strain is also reduced by a factor , due to cross- 
immunity. We assume that cross-immunity is weaker than within strain partial immunity, 
refiecting the fact that the antigenic distance between clusters must be larger than the cluster 
size. Therefore, we take a < a^, and 5 = dx — o" is a measure of the difference between the 
two competing strains, in terms of the immune response that they trigger. 

For a model with two (symmetrical) strains we have found that, at the level of mean field 
approximation, the reinfection threshold is the boundary between two different regimes. For 
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similar strains {6 ~ 0.02), the behaviour of the model below threshold is strain replacement, 
and the behaviour on and above threshold is strain coexistence. For distinct viral strains 
{6 ~ 0.2), the outcome is always strain coexistence, but the density oscillations are negligible 
above threshold, and very pronounced below threshold, implying that taking into account 
that the population is discrete, for realistic population sizes we find extinction of both strains 
below threshold ItII. This scenario has been confirmed by simulations of a stochastic model 

n 

for homogeneously mixed populations [16||. In the following we consider the contact structure 
of the host population and investigate its effect on pathogen diversity and evolution. 

We performed individual based simulations starting from an initial condition where the 
system is close to the steady state for a single resident strain and introduced a small frac- 
tion of individuals infected by the invading strain. We found that as local effects become 
important (low p), strain replacement (and also total extinction) are favoured with respect 
to the homogeneously mixed regime {p = 1). This happens both for similar {6 ~ 0.02) 
and for distinct strains {6 ~ 0.2), and, as opposed to the effect reported in the preceeding 
section, does not require p to be close to the epidemic persistence transition, only that the 
competing strains have moderate transmissibilities (see Figure 2). 

Our results differ from those reported in |17| for a model of strain evolution, with a 
high rate of supply of naive susceptibles, on a (static) small world network of contacts, 
where spatial correlations were found to favour pathogen diversity by reducing the spread of 
acquired immunity. In our model, with a low rate of supply of naive susceptibles, the main 
effect of spatial correlations is the screening of infectives that leads to a reduction of the 
effective transmissibility. As p decreases, this effect brings a system above the reinfection 
threshold closer to or even below the reinfection threshold, also favouring replacement for 
similar strains and either replacement or extinction for distinct ones. 

Since screening of infectives occurs quite generally in structured populations, the effect 
reported here should be fairly independent of the network's detailed topology and correla- 
tions. 
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IV. CONCLUSIONS 



Our results for the SEIR model support the conclusion that, for some purposes, successful 
modeling of disease spread must take into account that populations are finite and discrete, 
and must include a realistic representation of the spatial degrees of freedom or, more gen- 
erally, of the interaction network topology. We suggest a reassessment of the impact of 
environmental forcing by studying its effects on a realistic autonomous model as the one 
that we propose. In particular, fluctuation enhancement by spatial correlations might re- 
move some of the constraints on the strength of seasonality and avoid fine tuning of seasonal 
forcing amplitudes. 

In the context of multiple strain infections and the evolution and diversity of pathogens, 
the results of our simulations strongly support the conclusion that in systems with reinfec- 
tion and low rates of supply of naive susceptibles, the major effect of the host population 
structuring is a reduction of pathogen diversity. 
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FIGURES 

FIGURE 1 Infectives time series for homogeneously mixed and spatially structured 
populations from SEIR simulations on = 1000 x 1000 lattices for a) measles (tj = 8 days, 
Ti = 6 days) and b) pertussis (tj = 18 days, ti = 8 days), a) Grey line: Results for a 
homogeneously mixed population, p — 1, and transmission rate (3 — 2.4 day Black line: 
Results for a network with p — 0.2 and /3 = 4.75 day~^ also exhibiting an average period 
close to two years but much larger incidence oscillations, b) From top to bottom: Results for 
a homogeneously mixed population, p = 1, and transmission rate /3 = 1.5 day Results 
for a network with p = 0.2 and j3 = 4.0 day~^; Results for a network with p = 0.1 and 
(3 = 7.0 day~^. As for measles, the transmissibility is tuned to preserve the average period 
of the homogeneously mixed population. The increase in the incidence peak's amplitude 
is neghgible for p = 0.2, but for p = 0.1 the fluctuations are significantly enhanced, as for 
measles on a network with p = 0.2. 

FIGURE 2 Times series for invasion dynamics simulations of evolved viral strains, on 
small world networks over lattices with N — 800 x 800 nodes. The vertical axis measures 
the logarithm of the densities of infectives of the circulating and of the invading strains. We 

have taken 1/80 year^^ for the population birth rate, 203 year^^ for the transmissibility, 0.27 
for the reinfection parameter and 52 ycar^^ for the recovery rate of both diseases. Wc have 
considered, in (a), similar strains, with 5 — 0.02, and, in (b), antigenically distant strains 
with 6 — 0.2, and taken p — 1 (resp. p — 0.3, p — 0.1) in the first (resp. second, third) 
column. 
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